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NOMENCLATURE 


A 


ALPHA 


complex amplitude of the unsteady airfoil motion: 

for pitching, A * oscillatory angle of attack in radians; 
for plunging, A * displacement normalized by bslf-chord. 
The physical motion is Re(Ae^^-) 

mean angle of attack, deg 


C 


chord of airfoil, m 


CL 

CL, A 
CM 

CM, A 

CPU,A(CPL,A) 

CPU(CPL) 


exp(io)t) 
f ,FREQ 

IU,A(Q)(IL,A(Q)) 


mean lift coefficient, + up 

normalized unsteady lift coefficient, + up 

mean moment coefficient at leading edge, + nose up 

normalized unsteady moment coefficient at leading edge, 
nose up 

complex amplitude of the unsteady upper (lower) surface 

mean value of upper (lower) surface pressure coefficient, 
PU(?L) - PINF 
QINF 

cos o)t + i sin wt 
frequency, Hz, 

Qth moment of the complex amplitude of the unsteady upper 
(lower) surface pressure coefficient 


IU(Q)(IL(Q)) 

k,K 

M 

00 

P 

PINF 

PL,PU 

PTOT 

QINF 


Qth moment of the mean value of upper (lower) surface 
pressure coefficient 

reduced frequency, ^ 

free-stream Mach number 

complex amplitude of the unsteady pressur ; 
the physical pressure = Re(Pe^‘*^^) 

free-stream static pressure, N/m^ 

mean value of surface pressure, N/m^ 

total ’ essuie, N/m^ 

dynamic pressure, N/m^ 

PRECEDING \ 


v 



Re, RE 
T 
t 
U 

X,x 

a 


ot(t) 


chord Reynolds number 
period of the motion, sec 
time, sec 

free-stream velocity, m/sec 
distance along airfoil, m 

complex amplitude of unsteady angle of attack, deg 
mean angle of attack, deg 
instantaneous angle of attack, deg 


Complex notation: 
Im'' ] 

Mag[ ] 

Ph[ 1 

Re[ ] 


imaginary part of [ ] 
magnitude of [ ] 
phase of t ] > deg 
real part of [ ] 



EXPERIMENTAL UNSTEADY AERODYNAMICS OF CONVENTIONAL 


AND SUPERCRITICAL AIRFOILS 
Sanford S. Davis and Gerald N, Malcolm 
Ames Research Center 


SUMMARY 


Experimental data on the unsteady aerodynamics of oscillating airfoils 
in transonic flew are presented. Two O.S-m-'Chord airfoil models an 
NACA 64A010 and an NLR 7301 — were tested in the NASA-Ames 11- by 11-Foot 
Transonic Wind Tunnel at Mach numbers to 0.85, at chord Reynolds numbers to 
12 x 106 , and at mean angles of attack to 4®. The airfoils were subjected 
to both pitching and plunging motions at reduced frequencies to 0.3 (physical 
frequencies to 53 Hz) . 

The new hardware and the extensive use of computer-experiment integra- 
tion developed for this test are described. The geometrical configuration 
of the model and the test arrangement are described in detail. Mean and first 
harmonic data are presented in both tabular and graphical form to aid in com- 
parisons with other data and with numerical computations. 


1. INTRODUCTION 


The unsteady aerodynamics of both fixed- and rotary-wing airfoil sections 
must be thoroughly understood in order to provide safe margins for flutter, 
buffett, and other undesirable aerodynamic phenomena. This need is most 
apparent in the critical transonic speed regime where these detrimental 
effects are most prevalent. Recent developments in numerical simulations of 
transonic unsteady aerodynamics have also highlighted the need for new 
experimental activity in this area. In response to these needs, an extensive 
test program was developed at Ames Research Center to measure the unsteady 
aerodynamics of both a conventional and a supercritical airfoil under a wide 
range of flow conditions. The objective of the test was to measure unsteady 
pressure distributions at higher Reynolds numbers over a more extensive range 
of test conditions than had heretofore been attempted. This report presents, 
in graphical and tabular form, the mean and fundamental frequency data from 
that test. 

The data were obtained in the 11- by 11-Foot Transonic Wind Tunnel at 
Ames Research Center. Over 200 data sets, representing various combinations 
of airfoil geometry, Mach number, Reynolds number, mean angle of attack, 
motion mode, motion amplitude, and frequency arc reported. For each data 
set both the mean and first harmonic loads are tabulated, and the pressure 
distributions are presented in both tabular and graphical form. 



Section 2 describes the important features of the test apparatus in 
detail. Including the wind tunnel, model Installation, motion generators, 
model construction, and model geometry. (Some of the hardware was also 
described in ref. 1.) A discussion of the computerized data system, devel- 
oped especially for this test, is provided in section 3. The software was 
written such that on-line comparisons could be made between the current data 
set and theoretical predictions. The measuring system is also described in 
references 1 and 2. Section 4 outlines the test program and section 5 pre- 
sents the data. The method used to integrate the chordwise pressure distri- 
bution is described in appendix A, and the tabulated first harmonic pressure 
data, enclosed in microfiche form, is designated appendix B. 

Some of the data have already been analyzed and can be found in refer- 
ences 3-6. A small subset of the data has been selected by AGARD for inclu- 
sion in its ‘‘Standards for Aeroelastic Application it is cited in 
section 4. 


2 . TEST HARDWARE 


The arrangement of the apparatus and the special two-dimensional flow 
channel installed in the 11- by 11-Foot Transonic Wind Tunnel were based on 
the choice of an acceptable ratio of wind-tunnel height to wing chord 
(greater than 6). A chord of 0.5 m was chosen, resulting in the ratio 
(height) /(chord) = 6.8. Lowest hardware cost and minimum overall tunnel 
blockage could be obtained with a model spanning the tunnel, but construction 
of a full-span 0.5-m-chord model was impractical because first priority was 
assigned to obtaining high frequencies with minimal aeroelastic effects. An 
acceptable span-to-chord ratio of approximately 3 (1.33-m span) dictated the 
use of the splitter-plate arrangement shown in figure 1. Although previous 
investigators have successfully used splitter plates, a pilot test of the 
concept was nonetheless conducted in the Ames 2- by 2-Foot Tranoonic Wind 
Tunnel (ref. 7). T is test demonstrated that good quality transonic flow 
could be obtained with the chosen splitter plate arrangement - 

Figure 1 shows the general arrangement of the wing/splitter-plate/ 
actuator system as installed in the wind-tunnel test section. The normal 
3.35 m X 3.35 m test section was segmented with two steel splitter plates, 
3.35 m high by 2.8 m long. To minimize blockage, the thickness was the 
minimum necessary to accommodate the push-pull drive rods. To prevent exces- 
sive deflections of the splitter plates, side struts were instalJed for lat- 
eral support. The splitters extended into the tunnel's plenum area at the 
top and bottom; there they were bolted to I-beam anchors. Access panels for 
instrumentation cables aud clearance for the push-pull rods were included in 
the splitter plate design. 

The wing model was instrumented near its midspan station and attached to 
independently controlled hydraulic actuators through the push-pull rods. 

Thus, the wing was free to pitch and plunge in response to the actuator’s 
command signal. The wing was restrained in the fore-aft direction by a pair 
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cf carbon-epoxy drag rods, and in the lateral, roll, and yaw directions by 
sliding cover plates, which moved with the wing on the inner surface of the 
splitter plates. The hydraulic actuators, located in the lower plenum area, 
were supported by flexures; they bore directly into a massive concrete foun- 
dation through the four support columns. With this design, the tunnel pres- 
sure shell does not have to support the oscillatory reaction loads induced 
by the actuator's motion. 

The capabilities of the test apparatus include sinusoidal pitching 
oscillations over a frequency range of 0 to 60 Hz, with the maximum oscilla- 
tion varying from ±2° at low frequencies to ±0.8^ at 60 Hz about any chord- 
wise axis, and a vertical plunging motion up to ±5 cm (2 in.). 

The various components that make up the system will be described in more 
detail since the basic performance requirements dictated state-of-the-art 
designs in many cases. Many of the components are shown in the installation 
photograph in figure 2 and the pre-test setup in figure 3. In the following 
description it may be helpful to refer to these photographs to visualize the 
interrelationship among the various components. 


11- by 11-Foot Transonic Wind Tunnel 

The 11- by 11-Foot Transonic Wind Tunnel is a closed-return, variable 
density facility with a3.35x3.35x6.7m (11 x 11 x 22 ft) test section 
enclosed in a 6-m (20-ft) diameter cylindrical pressure cell. The air is 
driven by a three-stage, axial-flow compressor powered by four induction 
motors with a maximum continuous combined output of 135 MW (180,000 hp) . 

The Mach number can be varied continuously from 0.4 to 1.4 with the stagna- 
tion pressure variable from 50 kN/m^ to 225 kN/m^ (0.5 to 2.25 atm) resulting 
in Reynolds numbers rrom 6x10^ /m to 31x10^ /m. Maximum Mach and Reynolds 
numbers for this test were 0.85 and 25x10^ /m, respectively. 

The ventilated wall of the 11-Foot Transonic Wind Tunnel has a baffled 
slot arrangement (fig. 4). Six slots — 1.78 cm (0.7 in.) wide — between the 
splitter plates yield an effective open area ratio of approximately 8%. A 
resistive baffle fabricated from 0.16 cm (1/16 in.) sheet stock was inserted 
in each slot. The baffle is flush with the floor and ceiling, extends 
5.72 cm (2.25 in.) into the slot, and has a "wavelength** of 3.43 cm (1.35 in.). 


Splitter Plates 

Vertical splitter plates with trailing-edge flaps and horizontal side 
struts form the support structure for the wing. They each have a sharp lead- 
ing edge and a movable trailing-edge flap which is manually adjustable between 
±2® from the plane of the splitter plate. All testing was done with the flaps 
set at 0*^. Horizontal side struts attach to the outside of the splitter 
plates just below the horizontal plane of symmetry and protrude through the 
test section into the exterior structure. They provide stabilization and 
eliminate excessive lateral deflection from the aerodynamic loads. The 
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splitter plates were installed with a 0.1® divergence angle from tunnel 
centerline to account for boundary-layer growth. The thickness of the 
splitter plates varies in the streamwise direction in the following manner: 
following the sharp leading edge the next immediate section is 3.2 cm 
(1.25 in.) thi^k; it is followed by a 5-cm (2-in.) thick section in the 
center to accommodate the push-pull rods. The trailing-edge section is 
4.4 cm (1.75 in.) thick and tapers to a sharp trailing edge. The inside 
surface of the splitter plate is straight with all thickness variations tak- 
ing place on the outer surface. 

Openings in the splitter plate (figs. 5, 6) permit the wing to be 
attached to the top of the push-pull rods, which are centered in four channels 
cut into the lower portion of the splitter plates. When the wing is oscil- 
lating, sliding covers (figs. 7, 8) attached to the wing seal the openings. 
The covers are made of graphite epoxy to reduce weight and are Teflon-lined 
for free sliding. 


The splitter plates contain a total of 125 static-pressure orifices 
distributed over the inside and outside surfaces of both plates. The inside 
orifices were utilized to select the proper channel Mach number and, in con- 
junction with the outer taps, were used to monitor the loading on the 
splitter plates. While testing, accelerometers on the trailing-edge flaps 
were used to sense any large or potentially destructive flutter motions such 
as might be produced from the oscillating flow behind the wing or naturally 
induced from the channel air flow. 


Wings and Push-Pull Rods 

Model geometry- Two airfoil sections were chosen for this test program — 
one a conventional airfoil (an NACA 64A010) and the other a supercritical air- 
foil (the NLR 7301). The two wing models — span 1.35 m (53.2 in.), chord 
0.5 m (19.685 in.) — were designed to withstand accelerations of 
2.3x10^ m/sec^ (230 g) and aerodynamic loads of 44,000 N (10,000 lb). Both 
airfoils were subsequently chosen for inclusion in the AGARD standard series 
of test cases for aeroelastic applications (refs. 8, 9). Photographs of the 
models installed in the wind tunnel are presented in figures 7 and 8. Due to 
expansion of the molds in fabricating the models, the actual airfoil sections 
were slightly thicker than their theoretical counterparts. To expedite 
numerical simulations, three sets of ordinates are presented — the measured 
ordinates, smoothed versions of th»^ measured ordinates from Olsen’s computa- 
tions (ref. 8), and the theoretical ordinates. Because the measured ordinates 
contain large variations in the higher derivatives that adversely affected 
some trial solutions. It is recommended that either the smoothed or theoreti- 
cal ordinates be used for computing. Computations using the theoretical 
ordinates were satisfactory for the flow conditions attempted. 

The measured and theoretical airfoil sections are shown in figure 9. In 
each case he measurements correspond to the thicker section. Data for the 
NACA 64A010 and NLR 7301 airfoils are presented in tables 1 and 2, respectively. 
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Model instrumentation- The wing is instrumented with static pressure ^aps 
and dynamic pressure transducers, all of which are located at approximately 
midspan. The dynamic pressure transducers communicate to the wing surface 
via a small orifice with a small volume cavity. Locations of the static and 
dynamic orifices in both wings are shown in tables 3 and 4. It should be 
noted that dynamic transducers were not installed in the lower surface of the 
NLR 7301 airfoil. The lower surface unsteady pressures were sacrificed on 
that airfoil for the sake of increased resolution on the upper surface. 

Static pressure tubes are routed from the end of the wing through a cavity 
in the splitter plate to the tunnel plenum chamber below, and out an access 
port to scanivalve-transducer units exterior to the tunnel shell. Dynamic 
transducers are mounted in the wing by inserting the transducer (2.36 mm 
diameter) in the end of a long plastic sleeve, which is, in turn, inserted 
into a cylindrical channel molded into the interior of the wing. The sleeve 
terminates at the center of the wing at the orifice communicating to the 
wing surface. The lead wires are then routed out the opposite end of the 
sleeve in the wing (fig. 6) through the splitter plates and out through the 
tunnel walls to the data acquisition equipment in the tunnel control room. 

A single reference pressure tube from each dynamic transducer is also inserted 
into the plastic sleeve and routed through the splitter plate to the scanivalve- 
transducer assembly outside the tunnel. T^e transducer reference pressure can 
be selected to be the static pressure of the adjacent static orifice on the 
wing or any other selected pressure (such as the tunnel static pressure) . Six 
accelerometers were mounted inside the wing, one at each of the attachment 
points of the four push-pull rods near the corners of the wing and two at mid- 
span near the leading and trailing edges. The actual motion of the wing can 
be determined from the accelerometer output and compared with the output of 
tie motion transducers located in the actuator piston rods. These data 
Si.owed fiat the wing motions were faithfully recorded by the motion transducers. 

Model support s 3 ^stem- The wing model, mounted between the splitter plates, 
is connected to the push-pull rods through ci^ecial flexure bearings. The push- 
pull rods are, in turn, screwed directly into the actuator pistons. Both the 
wing and the push-pull rods are fabricated from a lightweight graphite-epoxy 
material. A short discussion of the fabrication of the rods and wings is 
given later in this section. The push-pull rods, 0.0412 m (1.625 in.) in 
diameter, are each capable of a 22,000 N (5,000 lb) tension load. The flex- 
ures located between the push-pull rods and the wing are also designed for a 
22,000 N (5,000 lb) load. A pair of graphite-epoxy rods mounted to the wing 
with a flexure support and attached to the splitter plates forward of the 
wing provide a means of counteracting the drag loads (see fig. 5) ; each rod 
can withstand 6,700 N (1,500 lb). 

Model fabrication- The fabrication of both the wing models and the push- 
pull rods required an extensive development effort by the Ames Model Develop- 
ment Branch. The requirements for maximum strength, stiffness, and light 
weight suggested the use of composite fiber materials. The problem of con- 
structing the wing was compounded by the requirement for internal mounting of 
the pressure transducers. The following description will illustrate briefly 
the steps used to fabricate the wing models. 
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The outside contour of the wing was defined by using a steel female mold 
split into two halves and machined to the coordinates of the airfoil section. 
The mold was made of common steel with a coefficient of thermal expansion 
compatible with that of the carbon graphite composite fibers. Using this 
mold, a fiberglass mold was made to construct the four interior silicone 
mandrels about which the carbon graphite fibers were wrapped to fcrm the 
internal structure. These mandrels were later removed from the wing after 
curing, leaving hollow sections between the webs (see fig. 6), The core 
mandrels were pregrooved for instrumentation tubes before wrapping with the 
appropriate number of layers of fiber strips. Tapered steel rods were inserted 
in the pre-cut grooves. The ribs between the four cores joining the upper and 
lower surfaces together were individually laminated and placed between the 
wrapped cores. A fixture was built to hold the four cores in place for wrap- 
ping of the entire model. Thirty-two plies of unidirectional graphite tape 
were wrapped around the pre-wrapped cores with plies at 0*^, 45®, and 90® with 
respuct to the chord. The model layup was then sandwiched between the two 
halves of the steel mold and caps were bolted on the ends. The entire assem- 
bly was heated in an autoclave o 250® F to expand the silicone core. This 
cycle forced the layup tightly against the interior walls of the mold. The 
model was subjected to a cure cycle of 350® F for 2 hr. After cooling, the 
model was removed from the mold and the cores removed. The tapered steel 
rods were also removed and a 0.5-mm hole drilled through the surface of the 
model to intersect the cavity left by the rods. The dynamic transducers were 
mounted in the end of a long plastic tube, which was inserted into the hollow 
cavity in the wing. The transducer body sealed against a shoulder in the tube 
forming a pneumatic seal. The volume between the transducer diaphragm and 
the orifice on the model surface was quite small, thus providing good high- 
frequency response. 

The push-pull rods were constructed of carbon-graphite fibers, using 
similar procedures. A two-part mold was made from mild steel with outside 
dimensions of 10.16 cm (4 in.) square by 203 cm (80 in.) long and a 4.13-cm 
(1.625-iii.) bore. A silicone core approximately 2.5 cm (1 in/^ in diamecer 
was made, and graphite fibers were carefully wrapped around it to a wall 
thickness of approximately 0.63 cm (0.25 in.). After curing and removing the 
core, the ends of the rocs were attached to steel end caps which provided an 
attachment point for the actuator system and flexures. 


Motion Generators 

The servo-hydraulic actuator system was designed especially for this 
test. It is driven by two 11-kW (150-hp) hydraulic pump units rated at 
4.1 X 10"3 mVsec (65 gal/min) at 20.7 x 10^ N/m^ (3,000 Ib/in.^). Each of 
the four actuators consists of two separate pistons on a single rod enclosed 
in a dual-chamber cylinder. The upper piston is used for generating dynamic 
forces, the lower piston for load biasing. The load bias system is necessary 
to support the mean aerodynamic lift load, thereby reducing the power required 
to drive the dynamic piston. As static bias requirements change, the servo- 
valve system responds to produce the required force output to maintain the 
set position. Velocity and fv>-ition transducers are combined, into a single 


6 



physical unit v^ith coils and cores aligned axially for mounting in the center 
of the actuator. 


Pretest Verification of System Components 

Because every part of this system was new, there was no test information 
available for judging performance and reliability of the apparatus. There- 
fore, a special pretest facility was built to permit a detailed checkout pro- 
gram. Many of the components, including the wing, push-pull rods, drag 
restraints, and the hydraulic actuator motion generator system, were new 
designs and could net be risked in the wind tunnel without pretest experi- 
ments. A special test stand v^as built for system verification. Figure 3 is 
a photograph of the assembly in the test area. A support structure was con- 
structed to which the various components were attached. The hydraulic actua- 
tors were mounted at the base with the push-pull rods attached to the top of 
the pistons. The wing was mounted on the push-pull rods with flexures and 
angle-of-attack blocks between the rod end and the wing end cap. The drag 
restraint was fastened on top of the rear flexures and the other end tied to 
the support frame. Lift loads were simulated by an inflatable bag between 
the lower surface of the wing and a support cradle fastened to the support 
stand. Drag loads were simulated by a pneumatically activated piston coupled 
to cables and straps looped over the wing. A nearly complete envelope of 
test conditions could be evaluated on the test stand. In the early stages 
cf the test checkout a wing constructed of fiberglass (shown in fig. 3) was 
used before risking the graphite-epoxy test wing. This proved to be an 
extremely valuable and low-risk method of evaluating the performance of the 
entire system. The only real limitations were that the fiberglass wing was 
not stiff enough to prevent large deflections at midspan (particularly in 
plunging) at frequencies above 30 Hz, and was not strong enough to accept the 
maximum lift loads. A limited amount of testing was done with the carbon- 
epoxy wing before installation in the wind tunnel. 


3. DATA ACQUISITION SYSTEM 


In the past, multichannel unsteady aerodynamic data were acquired using 
analog tape recorders where raw data were recorded and stored for future 
analysis. On-line analysis was restricted to a few selected channels, using 
special instrumentation to extract usable data from the great mass of incom- 
ing data. These systems suffered from long time lags between acquisition and 
analysis and the probability of unknowingly recording spurious data. In the 
present test a new computational data acquisition and analysis system was 
developed for on-line display of steady and unsteady aerodynamic data. Fig- 
ure 10 depicts the main elements of the new system. It can graphically dis- 
play the first-harmonic pressure distribution (both magnitude and phase) due 
to arbitrary pitch-plunge motions of the airfoil along with the conventional 
static pressure distribution. At the user's option, an overlay of selected 
theoretical or experimental pressure distributions from computer-resident 
codes or from a dedicated data bank can be accessed. 
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The system comprises a Data General Eclipse Model S/200 minicomputer, 
a high-speed (500 kHz) multichannel analog-to-digital converter, a large 
capacity (92 Mbyte) storage device, and a graphics terminal. The software 
system consists of approximately 50 independent Fortran-coded programs. The 
independent programs are controlled by two executive programs: one for 

dynamic data, the other for static data. 


Dynamic Data Acquisition 

The raw data come from a variety of sensors, the two most important 
being the airfoil motion (the input) and the surface pressures on the model 
(the output). The same sinusoidal signal that drives the four-channel 
hydraulic actuator system, which in turn drives the four push-pull rods 
attached to the wing, is also used to trigger a pulse to initiate the unsteady 
data acquisition process. Once the actuator control system is adjusted to 
impart the desired motion to the wing, the motion of the four push-pull rods 
is continuously monitored and acquired along with the unsteady data. 

The dynamic signals from up to 41 miniature pressure transducers are 
amplified and filtered before they enter the analog-to-digital converter. 
Because the signal is periodic, it is possible to obtain good waveform sam- 
ples with minimum storage per data point by signal-averaging the data. Theo- 
retically, a periodic signal is completely defined by just one cycle of data 
(e.g., a 40-msec record is all that is necessary to characterize a 25-Hz 
periodic oscillation) . However, the experimental signal is usually so con- 
taminated by random pressure fluctuations due to wind-tunnel and model- 
induced turbulence that one s^ycle of data is not very useful. 

The signal-averaging technique is implemented as follows: the raw data 

are synchronized ':ith a pulse train which is triggered at the same phase 
position for each cycle of the airfoil’s motion. These timing relations are 
shown ' figure 11. At time the sample waveform is recorded for t sec- 

onds. f\t time tQ + nT the waveform is recorded again for t seconds. The 
process is repeated M times. These M samples, each being initiated by 
the phase-locked pulse, are then ensemble-averaged to obtain the averaged 
signal. In the current experiment t is chosen to be slightly greater than 
one period, n = 2, and M = 100 is sufficient for a good average. At the 
user’s option, the signal-averaged waveform and its Mth realization for any 
selected channel can be displayed on the graphics terminal. 

For on-line analysis, the first harmonic of the response is most useful. 

A simple Fourier analysis algorithm is implemented to extract the magnitude 
and phase information at the fundamental frequency. These data are displayed 
in tabular form on the grapM^-i unit within 30 sec of the termination of data 
acquisition. These data are usually sufficient to determine if the unsteady 
data acquisition process was successful, cf more on-lim inalysis is required, 
the f irst-harr onic data may be displayed graphi'^ally in pressure coefficient 
form. The magtiitude and phase of the chordwise pressure distributions on the 
upper and lower surfaces of the airfoil are displayed along with certain 
theoretical curves, such as (1) linear, incompressible small-disturbance theory 
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(Theodorsen function) and (2) linear » compressible small-disturbance theory 
(Possio integral equation solver). For time-efficient on-line analysis it 
does net seem feasible to include unsteady transonic codes on the current 
generation of minicomputers. 

Also available for comparison are the results of other Investigations 
(theoretical or experimental) which have been stored in the data bank. For 
comparing with NACA 64A010 data, the theoretical investigations of Magnus 
(ref. 10) are available. For the NLR 7301 wing, experimental data obtained 
at NLR- Amsterdam (ref. 11) are available. It is possible to obtain a com- 
parison between the current data and the selected theoretical-experimental 
overlay in approximately 45 sec after the termination of data acquisition. 


Static Data Acquisition 

The static pressures are sensed with a conventional system using pneu- 
matic tubing connected to a pressure scanning valve. The electrical output 
of the pressure cell to which the unknown pressures are multiplexed are read 
with a digital voltmeter whose BCD (Binary Coded Decimal) output feeds 
directly into the minicomputer. 

The splitter-plate arrangement used for the oscillatory airfoil test 
requires special attention with regard to the free-stream Mach number (W^) • 

As discussed in previous reports (refs. 2, 7) the Mach number in the channel 
between the plates is not the same as chat computed from a static pressure 
tap in the plenum chamber. To obtain the approach Mach number (M^o) , the 
splitter plates are equipped with 125 static pressure orifices distributed 
among 10 rows above and below the plane of the wing on the inner and outer 
walls of the splitter plates. These pressures are also sensed by the scanning 
system. The computed Mach numbers on the splitters are displayed on the 
graphics unit and the approach Mach number is selected interactively by fair- 
ing the graphics unit's horizontal cursor to the data. Using this procedure, 
the Mach number can be selected to an accuracy of ±0.002. 

Once the Mach number has been chosen» the static-pressure distribution 
on tne wing is displayed along with selected overlays. A static pressure 
distribution with overlays can be displayed in approximately 30 sec after the 
raw data have been acquired. 


4. TEST PROGRAM 


As mentioned in the introduction, the test program was designed to meet 
the following primary goals: (1) to expand the existing unsteady test envelope 

to higher Reynolds numbers, higher reduced frequencies > different modes, and 
more diverse mean flow conditions; (2) to overlap the existing data base 
wherever possible; and (3) to provide a data base for the computation of 
unsteady transonic flows. A wide range of static and dynamic parameters was 
investigated in meeting these goals. The selected parameters are listed in 
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table 5. All of the data presented here were measured without a boundary- 
layer trip. Of the thousands of possible combinations, a subset of approxi- 
mately 200 comprises the current test matrix. Each selected combination is 
identified by a unique dynamic index (DI). A complete list of the test pro- 
gram in ascending numerical order is presented in table 6. The data in sets 
of "frequency sweeps" according to airfoil type and motion are given in 
table 7 for the NACA 6AA010 airfoil and in table 8 for the NLR 7301 airfoil. 

In reference 9 a series of airfoils was designated as AGARD standards 
for validating computational methods. The two airfoil sections used in this 
experiment are included in the standard series. Certain preferred flow con- 
ditions were chosen for comparative purposes. Ten cases for the NACA 64A010 
are presented in table 8 of reference 9. The corresponding dynamic indices 
are listed below: 

Test Case 123456 7 89 10 

DI 7 29 51 52 53 55 57 49 65 12 

Table 12 of reference 9 gives 14 test cases for the NLR 7301. These test 
cases were selected from the data reported in reference 11 and do not corre- 
spond exactly with the current series. In particular, mean flow conditions 
at the supercritical design point were slightly different. In table 8 of this 
report, the NLR 7301 frequency sweeps designated by rows 1-8 are the experi- 
mentally determined shock-free design conditions for this airfoil in the Ames 
11-Foot Wind Tunnel. They should be used for assessing computational methods. 


5. DATA REDUCTION AND PRESENTATION 


The primary output data are the pressure distributions on the airfoil 
along with quantities derivable from them. The data reduction and scaling 
applied to the raw data are described in this section. The data include 
static pressure coefficients, integrated static loads, complex amplitudes of 
the first harmonic pressure coefficients, and Integrated first harmonic loads. 


Static-Pressure Coefficients 

The static pressure data were converted to coefficient form using the 
conventional scaling: 


CPU (PU - PINF)/QINF ) 

1 

CPU = (PL - PINF)/QINF | 

where PU and PL are the measured mean pressures (in newtons per square 
meter) on the airfoil, and both PINF and QINF are flow parameters. (All 
symbols are defined in the nomenclature list.) The static-pressure data are 
presented in tabular and graphical form along with the dynamic data to be 
discussed subsequently. 
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Integrated Static Pressures 

The chordwise static pressure data were integrated according to the 
following formulas: 


f. 


1 


IU(Q) * I (CPU(X/C)lt(X/C)^]d(X/C) 


IL(Q) - I (CPL(X/C)]t(X/C)'"Jd(X/C) 

CL = IL(0) - IU(0) 

CM = -[IL(D - 1U(1)1 


(2a) 


(2b) 


The I- integrals are the Qth moments of the static-pressure distributions. 
Note that the moment coefficient is defined at the leading edge, nose-up 
positive. Some difficulties were encountered with the usual trapezoidal- 
type numerical integration scheme. The intejjration method that was ultimately 
adopted is described in detail in appendix A. 


Each static pressure run is associated with a unique identif ication — the 
static index (SI). Table 9 associates a static index with each dynamic index. 
Table 10 presents the integrated upper surtace, lower surface, and total load 
on the airfoil for each static index. 


Dynamic Pressure Complex Amplitudes 

The dynamic pressure data needed some preliminary processing. The first 
step was to Fourier-analyze the time-history data up to its fundamental fre- 
quency component. The fundamental frequency component is interpreted as a 
complex number that indicates its magnitude and phase shift with respect to 
the input motion. Figure 12 shows the steps used in decomposing the time- 
history into its real and imaginary components. The complex amplitudes 
CPU,A(CPL,A) are the quantities presented in this report. The physically 
realizable first harmonic unsteady pressure time-history is given by 

CPUD = fMag(A)](Re(CPU,A)cos wt - Im(CPU,A)sin ujtl \ 

} (3) 

CPLD * (Mag(A) ](Re(CPL,A)cos wt - Im(CPL,A)sin ojt] | 

where the complex amplitudes of the pressure coefficients CPU,A(CPL,A) have 
been normalized by the amplitude of the input motion Mag(A) . The time- 
history of the input motion is (fig. 12) 

A = [Mag(A)Jcos u5t 

where A is interpreted as an angular quantity for pitching motion or a 
translational quantity for plunging motion. 
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The 209 sets of first haraooic data, arranged by dynaBic index (DI), are 
presented graphically (real and imaginary parts) in figure 13. The correspond- 
ing static pressure distribution is also shown for reference. The tabulated 
static and dynamic data are presented in the enclosed microfiche (appendix B) . 
Note that only upper surface dynamic data were measured on the supercritical 
airfoil (DI > 115). 


Integrated Dynamic Pressures 


The first harmonic data were integrated in the same miiriner as the static 
data (eq. (1)). 


lO.A(Q) 

1L,A(Q) 

CL, A 

CM, A 


r 

f 


CCPU.A(X/C) 1 [(X/C)^ld(X/C) 


[CPL,A(X/0 ] ( (X/C)^]d(X/C) 


IL»A(0) - IU,A(0) 
-tlL.A(l) - I0,A(1)J 


(4a) 


(4b) 


These complex numbers are converted to time histories in exactly the 
same manner as in equation (3). The sign convention, interpretation of the 
lift and leading-edge moment, as well as the integration scheme, are the 
same as used in the preceding subsection on integrated static pressures. 


The six quantities in equation (4b) are given in table 11 for the 
NACA 64A010 airfoj.1 and in table 12 for the MLR 7301 supercritical airfoil. 


6. SUMMARY OF RESULTS 


Unsteady pressure data from an oscillating airfoil experiment in the 
Ames 11- by 11-Foot Transonic Wind Tunnel were presented. The data covered 
a wide range of parameters including airfoil geometry, mean flow conditions, 
motion mode, and frequency. These experimental results should be useful both 
for validating new computational methods and as an aid in aeroelastic analysis. 

To aid in the interpretation of the data, detailed discussions were 
included on the tunnel installation, tunnel geometry, and airfoil contour. 

The novel model fabrication and the new experimental techniques that were 
developed especially for this test program were also discussed. 

The data, presented in tabular and graphical form, include measurements 
of the mean pressure coefficients and real and imaginary parts of the funda- 
mental (first harmonic) frequency unsteady pressure coefficients. The 
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I^r-^ssure coefficient data are also presented in integrated form to facilitate 
in..erpretation of parametric trends. 

The data show that the unsteady aerodynamic response can be sensitive to 
> 11 of the parameters considered in this experiment. For subcritical flows* 
the two most important parameters are Mach number and frequency. In the range 
of mild transonic flows* airfoil geometry is an additional important factor. 
Ft.nally* in the flow regime where strong shock-wave/boundary- layer interactions 
are important* Reynolds number becomes another important parameter. This pro- 
g ession into increasingly complicated flows is consistent with the theoretical 
IT ithods that are used to predict these flows. Linearized subsonic theory 
i icludes the effect of the flow parameters M and k* and transonic theory 
orrectly accounts for airfoil geometry. In the most complex flow regime* 
kavier-Stokes modeling will be necessary to correctly predict the unsteady 
viscous interactions. 
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APPENDIX A 


METHODS FOR INTEGRATING EXPERIMENTAL PRESSURE DISTRIBUTIONS 


The integration of a function that is defined at a discrete number of 
points is not a simple problem. If a smooth curve can be fit through one or 
more of the discrete points, the integration becomes simple. The difficult 
part is to choose the appropriate family of smooth curves. 

A simple example will best illustrate the problem. Consider a pressure 
distribution having the functional form 

rP(X) = X /yT 0 < X < 1 (Al) 

This is the distribution that would be computed from thin*airfoil theory. The 
area under the curve, defined as the integral of equation (Al) bet een the 
limits 0 to 1 is 6.283 (to four significant figures). 

Nov consider a routine application of the trapezoidal rule. (The family 
of curves is simply the straight line connecting successive data points.) It 
is most convenient to consider the trapezoidal rule with equally spaced incre- 
ments. A typical case is shown in figure 14, where the function is divided 
into 20 strips. The value at the leading edge was approximated by extrap- 
olating the slope at the first chord position backwards to the leading edge. 
The computed loads for a 20-strip integration was 5.546. Compared with the 
exact area, this represents an error of 11. 7X. 

In actual practice, the leading-edge singularity is ameliorated by 
leading-edge bluntness and the errors computed above are probable upper 
bounds. However, in oscillating airfoil experiments, these leading-edge 
suction peaks can be quite high. This problem was pointed out some time ago 
by Runyan et al. (ref. 12). 

It the problem were only one of square-root-type singularities, an ele- 
gant solution is available. Gaussian quadrature techniques have been devel- 
opea (ref. 13) to approximate definite integrals by the finite sum 

N 

W(X)f(X)dX= 5] W.f. (A2) 

i=l 

They have not been widely used because the value of the function must be com- 
puted at sample points that are usually Irrational numbers. Standard 

tables are available giving the sample points and weights Wj^ for a given 
weight function W(X) . One such method ~ the Gauss-Jacobi quadrature — has a 

weight function W(X) = yi - X /VT. Figure 15 shows the sample points needed 
for a 20-strip Gauss-Jacobi quadrature. The computed area is 6.283, essen- 
tially the exact value. Gauss-Jacobi quadratures have been used extensively 
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in a recent theoretical report on oscillating airfoils in wind tunnels 
(ref. 14). A serious defect in the quadrature method for transonic flows is 
clear from figure 15. Sample points are sparsely located in the central 
region of the airfoil. Transonic flows with discontinuous pressure distri- 
butions (shock waves) are not well approximated by a scheme with such large 
increments in the region of discontinuity. Numerical experiments with dis- 
continuous pressure distributions have confirmed that unacceptably large 
errors can result. 

The numerical integration scheme that was finally adopted was the one 
described by Woodward (ref. 15). This method rectifies the leading-edge 
singularity by a simple transformation of variables: 



The pressure distribution presented in figure 14 is shown in the primed 
coordinance system in figure 16. The curve is finite ever)^here and a simple 
trapezoidal rule with 20 intervals has approximately six elementary trapezoids 
in the first 10% of chord. (Compare with fig. 14 where only two intervals 
constitute the first 10% of the chord.) Higher resolution is not compromised 
by a coarse mesh width in the area of discontinuities. Extensive numerical 
experiments have confirmed the validity of this procedure. For example, 
computations with a 20-strip integration resulted in an area of 6.288. The 
accuracy of the integrated quantities has been augmented by performing both a 
20- and 40-strip integration and by using a Richardson's extrapolation 
(ref. 16) to increase the accuracy. 

Once the integration scheme is selected, the remaining computational 
problem is to choose an acceptable interpolation-extrapolation scheme to 
transform the physical pressure tap locations to the desired mesh points. 

The method adopted was a polynomial fitting method for interpolating between 
data points and a linear extrapolation method for predicting values very 
near the leading and trailing edges. 
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APPENDIX B 


TABULATED FUNDAMENTAL FREQUENCY DATA 


Refer to the enclosed microfiche (inside back cover) for the 209 sets of 
tabulated steady and unsteady pressure data* 
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TABLE 1.- GEOMETRY OF THE NACA 64A010 AIRFOIL SECTION 
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TABLE 4.- NLR 7301 AIRFOIL SECTION: ORIFICE LOCATIONS 
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TABLE 3-- RANGE OF PARAMETERS USED IN TEST PROGRAM 


j Paraiaeter 

Syi^ol 

Range of values 

» 

Static quantities j 

Airfoil geometry 


NACA 64A010, MLR >301 ] 

Free-stream Mach number 

M„ 

0.45. 0.50. 0.65, 0.70. 0.75, 0.80, 0.85 

Mean angle of attack, deg 


0. 0.37, 0.57, 2.5, 4 

Reynolds number 

Re 

2.5>«10^ to 12.6*10^, depending on 

Dynamic quantities 

Motion mode 


Pitching, plunging 

Pitching axis location 


0.25C, 0.40C. 0.50C 

Pitching amplitude, deg 

a 

±0.25, ±0.50, ±1, ±2 

Plunging amplitude, cm 

h 

♦1 

Reduced f requeue y 

k 

0.025. 0.05, O.IO. 0.15, 0.20. 0.25. 0.30 
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TABLE 6.- TEST PROGRAM ARRANGED IN ASCENDING ORDER 
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TABLE 6.- Concluded 
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TABLE 7.- TEST PROGRAM ARRANGED ACCORDING TO FREQUENCY SWEEPS: NACA 64A010 
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TABLE 8.- TEST PROGRAM ARRANGED ACCOR.LING TO FREQUENCY SWEEPS: NLR 7301 
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TABLE 9.- IDENTIFICATION OF STATIC DATA CORRESPONDING TO EACH DATA SET 
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TABLE 10.- STATIC PRESSURES: INTEGRATED VALUES 
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TABLE 11.- NACA 64A010 AIRFOIL: DYNAMIC PRESSURE DATA, INTEGRATED VALUES 
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Figure 1.** General arrangement of oscillating airfoil test apparatus In 
NASA-Ames 11- by 11-Foot Transonic Wind Tunnel. 
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Wing model and push-pull mechanism 
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Figure 4.- Ventilated wall geometry of NASA-Ames 11- by 11-Foot Transonic 

Wind Tunnel. 


47 












Figure 3,- Detail of drag restraint and side strut support. 


Figure 6.- Wing end section with dynamic instrumentation leads. 
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Figure 7,- NACA ditAOlO model installati 



Figure H.- Ni.R 7301 model installation. 
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Figure 9." Sketch of airfoil profiles used In test program. 
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DATA FLOW 
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Figure 10.- Block diagram of the data acquisition scheme for the oscillating 

airfoil experiment. 




MOTION, A AMPLITUDE 


H-t-H 

Figure Il«- Tittlng diagraa for dynaaic data acquisition. Upp^r trace: 
dynamic data signal, t, is slightly greater than one period. Lover 
trace: trigger for analog-to-dlgital conversion, T * period, n « 2. 
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Figure 12.- Deco«aposition of pressure time histories into first harmonic 
complex amplitudes, (a) actual motion, (b) fundamental frequency com- 
ponent, (c) real and Imaginary parts (amplitudes proportional to cosine 
and sine waves, respectively). 
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Figure 13.- Continued. 
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Figure 13 .- Coneinucd 
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Figure 13.- Continued 
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Figure 13.- Continued. 
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Figure 13.- Continued 








Figure 13.- Continued. 
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Fi£ ire 13.- Continued. 




Figure 13.- Continued 
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Figure 13.- Continued. 



Figure 13.- Conclnued 




Figure 13.- Continued 



Figure 13.- Continued. 
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Figure 13.- Continued 
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Figure 13.- Continued 






Figure 13.- Continued 
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Figure 13.- Continued 













Figure > Continued 
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Figure 13.- Continued. 
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Figure 13.- Continued. 
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Figure 13,- Continued 



Figure 13.- Continued 
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Figure 13.- Continued 
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Figure 13.- Continued. 








Figure 13.- Continued 




Figure 13 .- Continued. 






Figure 13.- Continued 




Figure 13.- Continued. 





Figure 13.- Continued. 
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figure 13.- Concluded. 
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Figure 16,- Numerical Integration using 
transformed variables and trapezoidal 
lie. 



